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We present a consistent treatment of heavy quarks for jet production in DIS at NLO 
accuracy. The method is based on the ACOT massive factorization scheme and dipole 
subtraction method for jets. The last had to be however extended in order to take into 
account initial state splittings with heavy quarks. We constructed relevant kinematics 
and dipole splitting functions together with their integrals. We partially implemented 
the method in a MC program and checked against the known inclusive result for charm 
structure function. 

1 Introduction 

There are two basic approaches to heavy quarks production in DIS. First is so called zero-mass 
variable flavor number scheme (ZM-VFNS), which treats heavy quarks as massless partons with 
corresponding parton distribution functions (PDF). This scheme is applicable when the hard 
scale (taken here as the virtuality of the exchanged boson Q 2 ) is much larger then the mass 
mQ of a given heavy quark Q. On the other hand, when Q 2 is of the order of mq, so called 
fixed flavor number scheme (FFNS) is applicable. It retains the full mass dependence in the 
coefficient function and there is no PDF for Q as to the leading power it cannot appear in the 
soft part. 

Increasing precision of the data forces us to control also the intermediate region of Q 2 . The 
methods that address this problem are called general mass schemes (GM) [TJ [21 [3l 0] . They are 
however formulated for inclusive processes only and similar method relevant for jets is highly 
desirable. 

In the following we briefly describe our solution to this problem [5j [6] . It is based on the 
ACOT massive factorization theorem [TJ [7] and massive dipole subtraction method (DSM) 
8, 9,[TU], which however had to be reformulated in order to match with the former. 

2 Dipole subtraction method with massive partons 

Consider NLO calculation of a cross section for producing n jets in lepton-hadron reaction. The 
LO cross section is schematically written as 
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where AT is a normalization factor, f a are PDFs, d$„ !a is n-parton phase space (PS) and Ai n ,a 
is a tree-level matrix element (ME) with n final state partons and one QCD initial state parton 
a. The jet function F na determines the actual observable and is realized by a suitable jet 
algorithm. It satisfies -F n +i,a = F n ^ a in the singular regions of PS. At NLO the corrections 
involve loop diagrams living on n-particle PS and additional real emission belonging to (n + 1)- 
particle PS. Both contain IR singularities which ultimately cancel, however the cancellation is 
non-trivial as the singularities have different origin. An elegant and exact solution is provided 
by DSM. One adds and subtracts an auxiliary contribution T> n a , such that it mimics all the 
singularities of A^ n +i, a and at the same time can be analytically integrated over singular regions 
of PS. To be more specific we have 
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where virtual corrections to |A4„ !a | 2 are symbolically denoted as A4n°a P ^ ■ The subspace leading 
to singularities is d<p a and fulfils PS factorization formula d(fr n +i. a — d& n ^ a £g) d(j) a . Thanks to 
the properties of T> n a and F n+ i a the first square bracket is integrable in four dimensions, 
while in the second, the poles resulting from integral J d(j) a T> n a are cancelled against the ones 
in Mn°a pS> ■ However, not all singularities cancel in this way - there are also collinear poles 
connected with initial state splittings of massless partons. They are removed by means of 
collinear subtraction term C n a . It has the form 



C n ,a = ^2F ab ® \M n ,i 



(3) 



where T a b are renormalized partonic PDFs, i.e. the densities of partons b inside parton a. 
For instance in the MS scheme Tab ( z ) — ~ \ f 4 Pab ( z ), where P a b (z) are standard splitting 
functions. 

Within DSM the dipole function is realized as a sum of contributions corresponding to single 
emissions with different combinations of "emitter" and "spectator" partonfj^J Each such term 
D has a general form 

D = Vc\M n ,a\\ (4) 

where V is so called dipole splitting matrix (in helicity space) and encodes the information 
about some of the singularities of A^ n +i,o- The matrix C corresponds to color operators for 
parton emissions, which act on the matrix element. The notation above is symbolic and means 
that both quantities are correlated in the color and spin space. For DIS, there are three different 
classes of dipoles D, depending on whether emitter or spectator are in the initial or final states. 
Here we are mainly interested in the case of initial state emitter and final state spectator as 
they contain factorization-related information. 

When heavy quarks are present, the above general picture remains the same. If however a 
massive parton takes part in a splitting process there is no collinear singularity. Nevertheless 
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there are IR sensitive logarithms which become harmful when the external scale becomes large. 
We shall refer to such terms as quasi- collinear singularities |12j and abbreviate as q-singularitics. 

The first step towards GM scheme for 
jets is to construct dipole functions control- 
ling q-singularities for initial state emissions. 
Moreover we want to have possibly massive 
initial states as it is allowed by the ACOT 
scheme. It was partially done in [5] (for 
Q Qs splitting), while in jTU] the splitting 
processes with heavy quarks are considered in 
the final states only. 

Let us look at the particular example. 
Consider the initial state g — > QQ splitting. 
Let us assign the momentum p a to the gluon, 
Pi to the emitted final state quark (or anti- 
quark), and pj to the spectator. Using these, 
we construct new momenta which enter A4 n ,a 
in (Q: pf = w (pf + pf) — unit becomes a 
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Figure 1: ACOT charm structure function cal- 
culated using MC implementation of our method 
(MassJet). The calculations are done for xb ~ 0.05 
and CTEQ5 PDFs. 



new final state and p^ = (w — 1) (p% + p 1 ^ 
(u— 1) p% becomes a new initial state. The variables u, w can be determined from on-shell 
conditions for pj and pai- Our dipole splitting function reads 
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where [i r is a mass scale needed in D = 4 — 2e dimensions. In this case V is just diagonal in 
helicity space. 

Consider next the integral of ^ over one-particle subspace. It can be convenietly expressed 
in terms of rescaled masses rjf = mf /2pj ■ p a for some parton I. In the limit of small tjq we get 
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We see that there is a term of the form P gq (u) log i~Iq which becomes harmful when the scale 
becomes large (in massless case there would be a pole 1/e). Similar terms appear also in other 
dipoles for the initial state emissions. 



3 General mass scheme for jets 

In the spirit of the ACOT scheme, the initial state q-singularities have to be factorized out. 
It is accomplished by C„. a term with partonic PDFs T ao calculated in a special way. Let us 
recall at this point, that the latter are defined as certain ME of light-cone operators and can 
be calculated order by order using special Feynman rules (see e.g. |13|). The results contain 
UV singularities which have to be renormalizcd, leading to evolution equations. For the present 
application we calculate T a b to one loop with full mass dependence and renormalize them using 
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the MS schem^] Since counterterms are mass independent in this scheme, we assure that 
hadronic PDFs evolve according to standard massless DGLAP equations. For instance we get 

,2 

(7) 



•FpQ (z) = ^T R (l-2z(l-z)) log-^-. 
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We have checked that the above procedure leads to IR safe dipoles in the limit of vanishing 
r/Q. Moreover, the results coincide with those Rcf. [10 in the MS scheme. 

In order to perform numerical tests we have partially implemented our method in a dedicated 
C++ program based on FOAM [13] . Using the program, we have calculated the charm structure 
function F 2 and compared it with semi-analytical calculation in the ACOT scheme. This 
exercise uses three dipoles and two collinear subtraction terms. The virtual corrections are 
taken from [TS]. We find that the soft poles are indeed cancelled by the corresponding poles 
coming from the integrated dipoles. Moreover, we find agreement with the semi-analytical 
calculation and observe that our result properly interpolates between the two limiting solutions 
of the ZM-VFNS and FFNS schemes, as depicted in Fig. [TJ 

Let us stress that the result is obtained by a numerical integration of a fully differential cross 
section, which provides a severe test on the implementation of our massive dipole formalism. 
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